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Abstract 

We investigate the linear stability of plane Poiseuille 
flow in 2D under slip boundary conditions. The slip 
s is defined by the tangential velocity at the wall in 
units of the maximal flow velocity. As it turns out, 
the critical Reynolds number depends smoothly on 
s but increases quite rapidly. 

1 Introduction 

No-slip boundary conditions are a convenient ideal- 
ization of the behavior of viscous fluids near walls. 
In real systems there is always a certain amount of 
slip which, however, is hard to detect experimen- 
tally because of the required space resolution. In 
high precision measurements Elrick and Emrich pj 
detected slip of the order 0.1% in laminar pipe flow 
with Reynolds numbers of 16 to 4300. The measur- 
ing error in was nearly as low as the fluctuations 
due to Brownian motion. Very recently, Archer et 
al. observed the existence of slip in plane lami- 
nar Couette flow with added polymers. 

We examine how the linear instability of the 
steady plane Poiseuille flow depends on the slip s 
defined by 
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where u wall is the tangential velocity at the wall and 
u max is the midstream velocity. As boundary condi- 
tions we adopt 
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where z is measured in units of the channel half- 
width. The slip s is implicitly determined by the 
parameter b > with s — > in the limit b — > oo. 

2 Orr-Sommerfeld equation with 
slip boundary conditions 

The continuity equation in two dimensions is most 
conveniently satisfied by introducing a stream func- 
tion ^(x, z, t) where x denotes the streamwise direc- 
tion and z the direction normal to the boundaries 
(see Fig. [l]). The velocity (u,w) is connected to ^ 
through 
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In terms of ^> the Navier-Stokes equations for plane 
Poiseuille flow in two dimensions read in dimension- 
less form 
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As usually, is decomposed in the stream func- 
tion ^6 of the steady profile and a Fourier ansatz in 
x-direction for the disturbance field with the wave 
number a: 



$(x,z,t) = ^ b (z) + 
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%(z,t). (5) 



However, with slip the basic flow is now given by 
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The linearized part 
Sommerfeld equation 



of (|) leads to the Orr- 
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Figure 1: Geometry of the basic flow with slip 
boundary conditions. 



3 Numerical method 

We determine the critical (neutral) curves in the 
parameter space of the Reynolds number R and the 
wave number a of the disturbance. The Reynolds 
number is based on the channel half-width and on 
the midstream velocity of the steady flow. 

The solution of the differential equation (0) leads 
to a generalized eigenvalue problem that we solve 
numerically as in |3j using up to 70 Chebyshev poly- 
nomials as basis functions. The critical curve is 
the set of points (R, a) for which the most critical 
eigenvalue has zero real part with all other modes 
decaying exponentially. 

4 Results 

In Fig. H we present the critical curves for differ- 
ent slips s. The critical Reynolds number R c is the 
lowest Reynolds number on the critical curve. We 



define also the slip s c by the corresponding normal- 
ized tangential velocity of the critical mode at the 
wall. The results are listed in Tab. |l|. 

Obviously, the critical Reynolds number depends 
continuously on s. However, there is, perhaps sur- 
prisingly, a strong increase both of R c and s c with 
increasing slip s. In the limit b — > oo, i.e. s — > 0, 
one gets the well-known value R c ~ 5772. 



s 


0% 


0.1% 


0.2% 


0.5% 


1% 


Sc 


0% 


0.9% 


1.8% 


4.5% 


8% 


R c 


5772 


5773 


5781 


5847 


6070 


S 


2% 


3% 


4% 


5% 


6% 


Sc 


21% 


31% 


39% 


47% 


55% 


Rc 


6960 


8600 


11060 


15310 


23230 



Table 1: Slip s c of the critical mode and criti- 
cal Reynolds number R c at different slips s of the 
steady flow. 
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Figure 2: Critical curves of plane Poiseuille flow for 
different slips s = 0%, 1%, 2%, 3%, 4%, 5% and 6%. 
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